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ABSTRACT 



Aims. It was shown in previous works the existence of weakly chaotic orbits in the plutino population that diffuse very slowly. These 
orbits correspond to long-term plutino escapers and then represent the plutinos that are escaping from the resonance at present. In this 
paper we perform numerical simulations in order to explore the dynamical evolution of plutinos recently escaped from the resonance. 
Methods. The numerical simulations were divided in two parts. In the first one we evolved 20, 000 test particles in the resonance in 
order to detect and select the long-term escapers. In the second one, we numerically integrate the selected escaped plutinos in order 
to study their dynamical post escaped behavior. 

Results. Our main results include the characterization of the routes of escape of plutinos and their evolution in the Centaur zone. 
We obtained a present rate of escape of plutinos between 1 and 10 every 10 years. The escaped plutinos have a mean lifetime in the 
Centaur zone of 108 Myr and their contribution to the Centaur population would be a fraction of less than 6% of the total Centaur 
population. In this way, escaped plutinos would be a secondary source of Centaurs. 

Key words, methods: numerical - solar system: Kuiper Belt 



1. Introduction 

In the last few years the number of observed transneptunian ob- 
jects has enormously grown thanks to the progresses in the as- 
tronomical observations. This fact has allow to define, in a more 
rigorous way, the different dynamical classes previously iden- 
tified in the first years of discoveries. The transneptunian re- 
gion (TNR) can be structured into 4 dynamical classes (Chiang 
et al. I2007I) . The Resonant Objects are those in mean motion 
resonance with Neptune, the Classical Objects are those non- 
resonant objects with semimajor axis a greater than ~42 AU and 
low eccentricity orbits, the Scattered Disk Objects (SDOs) with 
perihelion distances q > 30 AU and large eccentricities, and 
the Centaurs Objects. This last group has perihelion distances 
inside the orbit of Neptune and they are transitory objects de- 
scendants of the other 3 classes, mainly from the SDOs, recently 
dislodge from the transneptunian zone by planetary perturba- 
tions (Levison & Duncan T19971 Tiscareno & Malhotra 2003), 
Di Sisto & Brunini 2007 ). Centaurs are sometimes defined ac- 
cording to their aphelion distance or to their semimajor axis (a), 
as for example the nomenclature of Gladman et al 2008 that uses 
a < a^ (where a^ is Neptune's semimajor axis), being then ob- 
jects entirely in the giant planet zone. Nevertheless, it is gen- 
erally accepted that they are objects which enter the planetary 
region from the TNR, evolve in the giant planetary zone and a 
fraction of them enter the zone interior to the Jupiter's orbit be- 
coming Jupiter Family Comets (JFCs). 

The resonant transneptunian population most densely pop- 
ulated are the plutinos, that are trapped into the 2:3 mean mo- 
tion resonance with Neptune, being Pluto its most representative 
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member. Some of the plutinos cross Neptune's orbit and hence 
one might think they would be sensitive to strong perturbations 
during close encounters with that planet. However conjunctions 
occur nearplutinos's aphelions and then close encounters do not 
occur. This fact provides a stable configuration of the resonance, 
with the critical angle cr = 3 A - 2A^ - rn librating around 180°, 
where A and m are the mean longitude and the longitude of per- 
ihelion of the plutino, and An is the mean longitude of Neptune. 

Duncan et al. ( 1995) analyze the dynamical structure of the 
transneptunian region through numerical simulations. They nu- 
merically integrate thousands of test transneptunian bodies in or- 
der to study their dynamical behavior and to determine which re- 
gions may be potential sources of the Jupiter family comets that 
we see today. With regards to the Plutinos, Duncan et al. ( 1995) 
showed that the boundaries of this long-lived mean motion res- 
onance have a time scale for instability of the order of the age of 
the Solar System, and so the long-term erosion of the particles 
is rather gradual and therefore must be continuing at the current 
epoch. The existence of these currently unstable orbits may be 
related to the origin of the observed Jupiter family comets. 

Morbidelli ( 1 19971) studied the dynamical structure of the 2:3 
resonance with Neptune in order to analyze possible diffusive 
phenomena and their relation to the existence of long-term es- 
cape trajectories. He performed numerical simulations integrat- 
ing the evolution of 150 test particles initially covering the 2:3 
resonance with eccentricities up to e — 0.3 and inclinations less 
than i = 5°, for 4 byr. The author found regular orbits that never 
escaped from the resonance present only at moderate eccentric- 
ity and small amplitude of libration, and a strong chaotic region 
at large amplitude of libration, which is quickly depleted and 
seems to be generated for the interaction between the vig and 
kozai resonances. Moreover, Morbidelli d 19971) showed the ex- 
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istence of weakly chaotic orbits that diffuse very slowly and fi- 
nally dive into the strong chaotic region. They are then long-term 
escapers and then plutinos recently escaped from the resonance. 
The source zone of those particles produce Neptune- encounter- 
ing bodies at the current epoch of the Solar System and should 
be an active source of Centaurs and comets at present. He found 
that 10 % of the bodies initially in that weakly chaotic zone are 
delivered to Neptune in the last 10 9 years. 

The first attempt at including the gravitational influence of 
Pluto into numerical models of the dynamics of plutinos was 
developed by Yu & Tremaine (11999b . These authors suggested 
that the effect of Pluto, ejecting objects from the 2:3 resonance to 
Neptune-crossing orbits, may contribute to or even dominate the 
flux of JFCs. At the same time, Nesvorny et al. (2000) analyzed 
the effect of Pluto on the 2:3 resonant orbits. They found that 
Pluto produces a large excitation of the libration amplitudes in 
the 2:3 resonance. However they estimated that the flux rate that 
contributes to the flux of short period comets is about 1 % of the 
2:3 resonant population per 10 8 years which is about the same 
value as the flux obtained by Morbidelli ( 1997) without Pluto. 

In the early 2000s, Melita & Brunini (2000) developed a 
comparative study of mean motion resonances in the transnep- 
tunian region. These authors used the frequency-map-analysis 
method (Laskar [19931 1 in order to describe the dynamical struc- 
ture of the 2:3, 3:5 and 1:2 Neptune resonances, at 39.5, 42.3 
and 47.7 AU, respectively. Particularly, Melita & Brunini (2000) 
showed that the 2:3 resonance presents a very robust stable zone 
primarily at low inclinations, where a great number of the ob- 
served plutinos are distributed. They suggested that the existence 
of plutinos in very unstable regions can be explained by physical 
collisions or gravitational encounters with other plutinos. 

de Elia et al. (2008) performed a collisional evolution of 
plutinos and obtained a plutino removal by their collisional evo- 
lution of 2 plutinos with R > 1 km every 10000 years or a flux 
rate of escape of 0.5 % of plutinos in 10 10 years. 

Very recently, Tiscareno & Malhotra (2009) carried out 1- 
Gyr numerical integrations to study the characteristics of the 2:3 
and 1:2 mean motion resonances with Neptune. Their main re- 
sults include maps of resonance stability for a whole range of 
eccentricities and inclinations. They made integrations with and 
without Pluto, and concluded that it has only modest effects on 
the Plutino population. They calculated the fraction of remaining 
plutinos in the resonance as a function of time and extrapolated 
this fraction after 4 Gyr and also evaluated the fate of escaped 
particles. 

From those previous works, plutino removal by "dynamic" 
is much greater than plutino removal by collisions. We refer to 
"dynamical", those numerical simulations that take into account 
the gravitational forces to follow up the evolution of a particle 
and occasionally cause the removal. Then in this work, we per- 
form "dynamical" numerical simulations in order to describe and 
characterize the routes of escape of plutinos and their contribu- 
tion to the other minor bodies populations of the Solar System, 
specially to Centaurs. 

2. The numerical runs 

The goal of our work is to characterize the post escaped evo- 
lution of plutinos and the presence of plutinos in other small 
body population in the current Solar System. So we need to 
identify the plutinos that have recently escaped from the reso- 
nance. Morbidelli (1997) studied the dynamical structure of the 
2:3 mean motion resonance with Neptune in order to analyze 
possible diffusive phenomena and their relation to the existence 



of long-term escape trajectories from the 2:3 resonance. The au- 
thor showed the existence of weakly chaotic orbits that diffuse 
very slowly finally diving into a strong chaotic region. These or- 
bits corresponds then to long-term escapers i.e. plutinos recently 
escaped from the resonance. Then we divided the numerical sim- 
ulations in two parts. At first we develop a numerical simulation 
of plutinos in the resonance in order to detect those plutinos that 
have recently escaped from the resonance. Second we perform 
a numerical simulation of the selected escaped plutinos in order 
to study their dynamical post escape behavior. In the following 
subsections we will describe both integrations. 

2.1. Pre-runs. The integration in the resonance 

In order to detect the long-term escapers from the plutino popu- 
lation, we performed a numerical integration following the study 
of Morbidelli ( 1997). We integrate the evolution of 20, 000 test 
particles under the gravitational influence of the Sun and the four 
giant planets over 4.5 Gyr with an integration step of 0.5 years 
using the hybrid integrator EVORB (Fernandez et al. 2002). 

We set the initial orbital elements such that they cover the 
present observational range of orbital elements of plutinos. The 
initial semimajor axis of the particles was set equal to the ex- 
act value of the resonance a-, - 39.5 AU. The initial argument 
of perihelion a>, longitude of node Q and the mean anomaly M 
have been chosen at random in the range of [0°, 360°] in a way 
that the critical angle <x remains between 180° and 330°. Since <x 
librates around 180°, and given the relation between a and <x in 
the resonance, the election of <x > 180° and a, = 39.5 AU. covers 
the 2:3 mean motion resonance (see Morbidelli 1 19971 for a com- 
plete explanation). The maximum limit of cr equal to 330°, was 
selected taking into account previous papers of Malhotra ( 1996) 
, Morbidelli 09971 ) and Nesvorny and Roig (2000), that state 
that orbits starting at large amplitude of libration (A^ ) are in a 
strong chaotic region, very unstable and fast driven to the bor- 
ders of the resonance. In the present paper we are interested in 
the long term escapers from the plutino population and because 
particles with large initial A a will escape at the begining of the 
integration they will not contribute to the long term escaper flux. 

The initial eccentricity and inclination of the particles have 
been randomly chosen in the intervals [0,0.35] and [0°,45°], re- 
spectively. 

The test particles were integrated up to the first encounter 
within the Hill sphere of a giant planet, collision onto a planet 
or ejection. Those cut off conditions mean that plutinos moved 
from the stable zone of the resonance, so that they suffer an en- 
counter with a planet, and then they do not belong any more to 
the resonant population. Therefore they are the escape condi- 
tions of plutinos. 

2.2. Selection of the long-term escapers 

From the 20, 000 initial particles, 17, 577 (87.9%) left out of the 
integration at a given time either because of an encounter with a 
planet or an ejection. We have 21 particles that are ejected of the 
Solar System and 17, 556 that encounter Neptune or Uranus. The 
remaining 2,423 particles (12.1%) keep inside the resonance up 
to the end of the integration. We will call the particles that leave 
the resonance, either because of an encounter or ejection, "es- 
caped plutinos". From our simulation we can calculate the rate 
of escape of particles from the resonance. In Fig. Q] we plot the 
cumulative number of escape particles from the resonance (N e ) 
with respect to the number of the remaining particles N p , where 
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N„ 



20,000 - N e , as a function of time. It can be seen that 



the number of escaped particles raise quickly at the beginning 
up to t ~ 1 .5 Gyr. In this point, the slope of the curve changes 
and behaves roughly as a linear relation. This change of slope 
was already noticed by Morbidelli (1997) and is related with 
the time when the strongly chaotic region is completely depleted 
and the weakly chaotic region starts to be the dominant source 
of Neptune-encountering bodies. We fit to the plot for t > 1.5 
Gyr, a linear relation given by: 



N e /N p - at + b, 



(1) 



9.00713xl0- 10 +4.735xl0~ i y- 1 and/? : 



3.28079; 



where: a = 
0.001347. 

Morbidelli d 19971) sets out that the slow diffusion region is 
the only active source of the 2:3 mean motion resonance that 
produce Neptune encountering bodies at present; then the slope 
of the linear relation, a, represents the present rate of escape of 
the particles from the 2:3 mean motion resonance. We have no- 
ticed that there are very few particles of which the semimajor 
axis diverge from the resonant value before they have an en- 
counter with a planet or even if they never have an encounter 
with a planet. But this behaviour doesn't change the previous 
calculated rate of escape from the resonance. 

On the other hand, we plot in Fig. [2 the number of re- 
maining plutinos (N p ) as a function of time. It can be seen 
a change of slope a t t ~ 1 00 Myr, that was also already no- 
ticed by Morbidelli JI997) and Tiscareno & Malhotra (2009). 
We can fit a power law to the number of surviving plutinos 
versus time for t > 100 Myr and is given by: 



N„=k t f 



(2) 



where: k = 7560950 + 34600 and/3 = -0.362672 + 0.0002372. 

Previous works have used power laws for fitting the num- 
ber of surviving particles in the 2:3 mean motion resonance 
(Morbidelli 1997, Tiscareno & Malhotra 2009). In particular 
treating Morbidelli (1997) data (see his Fig. 12) for t > 100 
Myr one can obtain /? = -0.32. This is a value very close 
to our value of p - -0.36. As for the fitting of Tiscareno & 
Malhotra (2009), they obtained a somewhat steeper slope, al- 
though they fit the number of remaining particles vs time for 
the last 0.5 Gyr of their integration (this is from 500 Myrs to 
lGyr). 

The number of small objects in the plutino population is not 
well determined, since the present observational surveys can't 
cover all the small-sized objects. Then, the size distribution of 
plutinos is calculated from different surveys up to a given size, 
typically up to a radius R ~ 30 km. For objects with radius less 
than this size, the population could have a break (Kenyon et al. 



2008 Bernstein et al. 2004 Elliot et al. 2005). There are also the- 
oretical models which account for accretion and fragmentation 
of planetesimals in the region of the Kuiperbelt and also predict 
a broken power-law size distribution at a given radius (Kenyon 
etal. |2008). 

de Elia et al. ([2008 ), analyzed the size distribution of pluti- 
nos, taking into account the suggested mass of the population 
and the possible existence of a break in the distribution. Through 
their collisional evolution of plutinos they concluded that the ex- 
istence of a break in the plutino size distribution should be a pri- 
mordial feature. Then considering the analysis made by de Elia 
et al. (120081 ) and the three plutino size distributions proposed, 
we calculate the present number of plutinos with radius R > 1 
km as N p ~ 10 8 - 10 9 , depending on the existence of the break. 




2000 2500 
t[My] 



4500 



Fig. 1. Cumulative number of escaped particles (N e ) with respect 
to the number of the surviving particles N p , versus the time t in 
Myr. The fit to the plot for t > 1 .5 Gyr is also showed. 
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Fig. 2. Number of remaining plutinos (N p ) versus the time t in 
Myr. A power law fit to the plot for t > 100 Myrr is also showed. 



Then we will have a present rate of escape of 1 to 10 plutinos 
with R > 1 km every 10 years. 

We have 1 183 particles that escape from the resonance after 
t — 1 .5 Gyr and so they come from the slow diffusion region of 
the resonance. Since those particles would represent the present 
escaped plutinos, we identified their original orbital elements for 
our second integration. We will explain this in the following sub- 
section. 



2.3. The Post-escape integration 

We numerically integrated again from t — the 1183 particles 
that escape at t > 1.5 Gyr., in order now, to obtain their Post 
escape evolution. We integrated those particles with the same 
computing conditions than the previous integration. We used the 
EVORB code under the gravitational influence of the Sun and 
the four giant planets, the same step of integration of 0.5 years, 
but now we followed the integration for 10 Gyr. The particles 
were removed from the simulation when they either collide with 
a planet or the Sun, or when they reach a semimajor axis greater 
than 1,000 AU or they enter the region inside Jupiter's orbit 
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(r < 5.2) where the perturbations of the terrestrial planets are 
not negligible. We recorded the orbital elements of the particles 
and the planets every 10 4 years. In the following sections we are 
going to describe the results of this simulation. 



3. Escaped plutinos. General results 

The great majority of escaped plutinos have encounters with 
Neptune, so this planet is the main responsible for their post es- 
cape evolution. We register the encounters with the major plan- 
ets at a distance of less than 3 Hill's radii. We have 1,945,643 
encounters with Neptune, 534,557 encounters with Uranus, 
88, 947 encounters with Saturn and 1, 309 with Jupiter. However 
the encounters with Jupiter in particular are reduced by the re- 
moval of objects at a distance of Jupiter. 

From the 1183 initial particles, 1179 are removed from the 
integration and 4 particles remain in the integration. Those four 
particles have encounters with Uranus or/and Neptune, they have 
short incursions to the Centaur zone and afterwards they are 
quickly transferred to a mean motion resonance in the SD re- 
maining there up to the end of the integration. From the 1179 
particles removed from the integration at some time, 790 (67%) 
are ejected, 385 (32.7%) reach the zone interior to Jupiter's orbit 
(r < 5.2), we have 4 (0.3%) collisions with the planets: 1 with 
Saturn, 1 with Uranus and 2 with Neptune. Those numbers can 
be compared with the ones obtained in the numerical simulation 
performed by Di Sisto & Brunini (2007); in particular the num- 
ber of escaped plutinos that reach the orbit of Jupiter and the 
number of ejections are similar to those numbers for SDOs with 
low semimajor axis and perihelion distances less than 35 AU. 
Also we can compare our results with the ones by Tiscareno & 
Malhotra (2009). They obtained that 27% of escaped particles 
reach the zone of r < 5.2. This is a number slightly small than 
ours, but it could be due to the fact that Tiscareno & Malhotra 
(2009)'s result was obtained from the escaped particles from 
their 1 Gyr integration. 

When a plutino escape from the resonance, it is transferred 
to the Scattered Disk (SD) zone (q > 30 AU) or to the Centaur 
zone (q < 30 AU). This is also noticed by Tiscareno & Malhotra 
(2009). The mean scale of time in reaching the Centaur zone is 
670, 000 years and in reaching the SD zone is about 6 Myr. 

3. 1 . Distribution of orbital elements 

In Fig. OJl we plot the time-weighed distribution of escaped 
plutinos in the orbital element space. There it is represented the 
probability distribution of finding an escaped plutino in the or- 
bital element space. These plots assume time-invariability, so 
they don't represent the real case where plutinos are continu- 
ously leaving the resonance, passing through a certain zone out 
of the resonance and leaving the solar system, but they help to 
identify the densest and empties regions. As we can see, Neptune 
is the planet that mainly governs the post escape evolution of es- 
caped plutinos. This behavior can be seen in fig <£3j as the densest 
zone near Neptune's perihelion. The densest zone in the orbital 
element space of escaped plutinos corresponds to the ranges of 
30 < a < 100 AU and 5° < i < 40°. Also it can be seen several 
mean motion resonances densely populated, those corresponds 
to the blue lines denoted in the a vs i plot. 

As we mentioned in the previous section, escaped plutinos 
take up the Centaur zone and the SDO zone. And in general 
they switch each of those populations during its dynamical evo- 
lution until it ends like one of them. This is obviously due to 



the presence of Neptune that lead the dynamical evolution of 
escaped plutinos. In particular, it is notable that when escaped 
plutinos are transferred to the SD they are quickly locked into a 
mean motion resonance with Neptune. This behaviour is simi- 
lar to the behavior of SDOs analyzed by Fernandez et al. (2004) 
and Gallardo (2006). The most densely populated mean motion 
resonances are: 2:3N, 4:7N, 4: UN, 1:3N, 1:6U and 1:5N, where 
"N" means "with Neptune" and U means "with Uranus". 

The distribution of escaped plutinos in the Centaur zone 
is similar to that obtained from SDOs by Di Sisto & Brunini 
(120071) . We will analyze this contribution in the next section. 

4. Contribution to the Centaur population 

4. 1 . Mean lifetime 

From the 1 1 83 particles that escape from plutinos 1 179 particles 
enter to the Centaur zone. The four remaining escaped plutinos 
are ejected by the first encounter with a planet (Neptune in all 
these cases) within 1 Hill's radii. The escaped plutinos have a 
mean lifetime in the Centaur zone of Iq - 108 Myr. This is 
greater than the mean lifetime of Centaurs from the SD of 72 
Myr (Di Sisto & Brunini 2007). In Fig. © we plot the normal- 
ized fraction of escaped plutinos against intervals of lifetime in 
the Centaur zone. Also it is plotted, for a comparison, the distri- 
bution of lifetimes in the Centaur zone for objects from the SD. 
As we can see Centaurs from plutinos have greater lifetime than 
Centaurs from SDOs, and the great majority of lifetimes of pluti- 
nos in the Centaur zone is greater than 1 Myr. So escaped pluti- 
nos spend long time in the Centaur zone, from 1 Myr to 1000 
Myr. From our numerical simulation, we noticed that the mean 
lifetime of plutinos in the Centaur zone, as a function of the ini- 
tial inclination, has nearly the same behavior than Centaurs from 
SDOs (see Fig.5 in Di Sisto & Brunini 120071 1. but they show 
higher values in each bin of initial inclination. So the difference 
in the mean lifetime of plutino-Centaurs and SDO-Centaurs is 
not dependent on the initial inclinations. 

In Fig (0 we plot the lifetime of plutinos in the Centaur zone 
with perihelion distances q less than a given value and also the 
mean lifetime of SDOs in the Centaur zone, as a comparison. 
For building this plot, we count for each particle (i) the time that 
it spends with q less than a given value qo, {dtjiq < qo)). Also 
we count the number Af of particles that remain for a certain time 
with q < qo. Then the mean lifetime of the particles with q < qo 

wouldbe/ c (<4o) = ^4^. 

We can see a stronger dependence of Centaur lifetime on 
their perihelion distance for plutino-Centaurs than for SDO- 
Centaurs. However the greatest differences are found for large 
perihelion distances, and for perihelion distances less than ~ 15 
AU both populations have similar lifetimes. Then, escaped- 
plutinos live more time than SDOs in the greater-perihelion 
Centaur zone, causing a slower diffusion to the inner Solar 
System of escaped-plutino orbits than of SDOs orbits. We sug- 
gest that this fact is due to the difference in the relative encounter 
velocities with the major planets of SDOs and plutinos. In Fig. 
(O we show the distribution of the relative encounter veloci- 
ties of escaped plutinos (solid line) and SDOs (dotted line) with 
the four giant planets. It is clearly seen that SDOs have smaller 
encounter velocities than plutinos. The mode of the encounter 
velocities is equal to 3.22 km/s for plutinos and 1.86 km/s for 
SDOs, and the mean values are equal to 2.92 km/s for pluti- 
nos and 2.45 km/s for SDOs. Orbits with smaller relative en- 
counter velocities suffer bigger changes in velocity, and then also 
in the orbital energy. Successive encounters with the major plan- 
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Fig. 3. Time-weighed distribution of escaped plutinos in the or- 
bital element space. 



ets make the objects to evolved to the inner Solar System or to 
the transneptunian region up to ejection. So bigger changes in 
the orbits lead to accelerate those transitions. We conclude that's 
why SDOs have a mean lifetime smaller than that of escaped 
plutinos. 



1e+04 1e+05 1e+06 1e+07 1e+08 1e+09 
l c [years] 

Fig. 4. Normalized fraction of escaped plutinos against lifetime 
in the Centaur zone (filled line) and normalized fraction of SDOs 
that enter the Centaur zone against lifetime there (dashed line). 



100 




q[AU] 

Fig. 5. Mean lifetime of Plutinos (filled line) and of SDOs 
(dashed line) in the Centaur zone with perihelion distances less 
than a given value. 




Fig. 6. Normalized distribution of encounter velocities with the 
major planets of SDOs (dashed line) and escaped plutinos (solid 
line). 
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4.2. Number of plutino-Centaurs 

In order to calculate the number of escaped plutinos located at 
present in the Centaur population, we calculate the present rate 
of injection of escaped plutinos into the Centaur zone. As we 
have mentioned we consider that the long-term escapers of the 
plutino population, i.e., those that escape after t = 1.5 Gyr, rep- 
resent the present plutino espapers. In Fig, (jT) we plot the cumu- 
lative number of escaped plutinos injected into the Centaur zone 
(N c ) with respect to the number of the remaining plutinos (N p ), 
as a function of time. We have taken into account the fact that 
a plutino escape from the population when it has an encounter 
with a planet or it is ejected, so the number of the remaining 
plutinos depends on plutinos injected to the Centaur zone but 
also on the escaped plutinos that are permanently or transitory 
injected into the SD zone. As Fig. (0 shows, the ratio N c /N p is 
well fitted by a linear relation, given by: 



NdN p - ct + d, 



(3) 



where: c = 1.62076 x 1(T 10 + 8.502 x lO^'V -1 and d = 
-0.228442 + 0.000242. 

The slope of the linear relation, c, represents the present rate 
of injection of plutinos from the 2:3 mean motion resonance into 
the Centaur zone. 

The present number of escaped plutinos in the Centaur pop- 
ulation can be calculated by differentiating Eq. 2: dNc/dt = cN p 
where we have assumed that N p is constant during an interval 
of time dt. Then, in order to calculate the present number of 
plutino-Centaurs we take dt - Iq and now assume that N p is 
constant during the lifetime of plutino-Centaurs. 



N c =cN p l c 



(4) 



where Iq is the mean lifetime in the Centaur zone. 

Taking from de Elia et al. (2008) the present number of pluti- 
nos with radius R > 1 km as N p ~ 10 8 - 10 9 ; the rate of injection 
of plutinos with radius greater than 1 km to the Centaur zone will 
be between 1.6 to 16 plutinos every 100 years. This is between 
250 and 25 times less than the rate of injection of SDOs to the 
Centaur zone obtained by Di Sisto & Brunini (2007). Then the 
present number of plutino-Centaurs with radius greater of 1 km 
would be between 1.8 x lO 6 - 1.8 x 10 7 . 

Di Sisto & Brunini (2007) estimated a number of Centaurs 
with R > 1 km coming from the SD of ~ 2.8 x 10 8 , then Centaurs 
coming from plutinos would represent a fraction of less than 6% 
of the total Centaur population. That is to say that the plutino 
population is a secondary source of Centaurs, comparable to 
the contribution of the low eccentricity transneptunian objects 
according to the estimations of Levison & Duncan ( 1997j> of 
1.2 xlO 7 . 



4.3. Orbital evolution of escaped plutinos in the Centaur zone 



Di Sisto & Brunini (2007) analyzed the dynamical behavior of 
the SDOs that enter the Centaur zone, and found four classes of 
dynamical evolution. The Centaurs, then, behaves as one of the 
classes or as a combination of them. For the reason of complete- 
ness, we will briefly describe the four classes here. The first type 
is characterized by the conservation of the perihelion distance 
in a range of values between that of Saturn and Neptune's orbit, 
the conservation or the very slow variation of the perihelion lon- 
gitude, and eccentricities greater than ~ 0.8. The characteristics 
of this evolution make the orbit go into a "pseudo-stable" state 
during which the encounters with the planets are avoided or very 
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Fig. 7. Cumulative number of escaped plutinos injected into the 
Centaur zone (N c ) with respect to the number of the remaining 
plutinos N p , as a function of time. The dashed line is the fitting 
to the data. 



weak, causing a very slow variation of the orbit orientation. The 
second type of objects are those that show "resonance hopping", 
(this is a phenomenon in which objects move quickly from one 
resonance, in this case with Neptune, to another) combined with 
a behavior similar to the first one but with less constant eccen- 
tricity values and constant perihelion distances for shorter inter- 
vals of time (Tiscareno & Malhotra. l2003l l. This type of objects 
have also transfers between mean motion resonances and Kozai 
resonances. A Kozai resonance occurs when the argument of 
pericentre, w, librates about a constant value. For low inclina- 
tions it is possible for w to librate about w — 0° and w = 180°, 
and for large inclinations about w = 90° and w - 270°. The 
semimajor axis of the object remains constant but the eccentric- 
ity and the inclination of the orbit are coupled in such a way that 
e is a maximum when i is a minimum, and vice versa. In these 
two first types, the objects have casual encounters with Neptune 
and Uranus, sometimes also with Saturn, but they are not strong 
enough to drastically change their orbit causing a kind of sta- 
ble orbit in the Centaur zone. In the third type of objects, we 
group those that have the behaviors of the first and second types, 
but they have perihelion distances near Neptune. So, the objects 
are continuously entering and leaving the Centaur zone. The last 
type of objects are those that enter a mean motion resonance or 
Kozai resonance for almost all their lifetime as a Centaur. 

Bailey & Malhotra (2009), have analyzed the chaotic behav- 
ior of the known Centaurs. Their analysis has revealed that two 
types of chaotic evolution are quantitatively distinguishable. One 
random walk-type behavior and an orbital evolution dominated 
by intermittent resonance sticking. These two dynamical classes 
embrace the four ones already found by Di Sisto & Brunini 
(2007). Bailey & Malhotra (2009) also found that these two 
types of behavior are correlated with Centaur dynamical life- 
time. 

In this paper we have analyzed the orbital evolution of es- 
caped plutinos in the Centaur zone, and we have found that they 
can be grouped into the four dynamical classes proposed by Di 
Sisto & Brunini (120071 ). There are more particles that have the 
dynamical behaviors of the second class. It is notable the great 
frequency of the presence of Kozai resonances in all the four 
classes. In the Centaur zone, mean motion resonances and Kozai 
resonances are more frequent for semimajor axis between 30 
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Fig. 8. Dynamical Evolution of the orbital elements of an es- 
caped plutino injected into the Centaur zone. 



AU and 50 AU. For example in Fig. ([8} we show the dynam- 
ical evolution of one plutino-Centaur of the second class that 
shows transfers between mean motion resonances and Kozai res- 
onances great part of the time. There are also some particles that 
once they escape, they return to a plutino state for some time and 
then continue their evolutions until their final state. Also there 
are particles that exhibit a "hand off" from the gravitational con- 
trol of one planet to another until crossing the Jupiter's orbit. 
The time scale of this hand off is whether as short as few Myr 
or as long as some hundreds of Myr. The particles that have the 
shorter lifetimes (say some Myr) carry out short incursions to the 
Centaur zone, or a handing down from a Jovian planet to the next 
inside until Jupiter, eventually passing through mean motion res- 
onances, or a quick diffusion through all the inter giant-planetary 
zone until ejection or injection to Jupiter orbit. 



5. Conclusions 

We have performed two numerical simulations in order to first 
obtain particles representative to the plutinos that are escaping 
at present from the resonance and second in order to describe 
their dynamical post escape evolution. In the first simulation we 
integrate 20,000 initial particles in the 2:3 resonance and find 
that ~ 88% of the particles left out the integration and the rest 
remain in it. Considering a plutino population with radius grater 
than 1 km of N p ~ 10 8 — 10 9 , we obtained a present rate of 
escape of plutinos between 1 and 10 every 10 years. From this 
integration we selected those particles that are representatives of 
the present escape plutinos and performed a second integration. 
From this last integration we obtained the dynamical evolution 
of plutinos once they escape from the resonance. From the 1 183 
initial particles, 1179 were removed from the integration and 4 



remain in it. From the 1 179 removed particles, 787 are ejected, 
385 reached the Jupiter's zone and 4 collide with the planets. 

We found that the great majority of escaped plutinos have 
encounters with Neptune, and this planet governs their dynam- 
ical evolution. When a plutino escape from the resonance, it is 
transferred to the SD zone (q > 30 AU) or to the Centaur zone 
(q < 30 AU) but it eventually switches to those population, due 
to the dynamical influence of Neptune. 

The densest zone in the orbital element space of escaped 
plutinos corresponds to the ranges 30 < a < 100 AU and 5° < 
i < 40° and perihelions near the orbit of Neptune. When escaped 
plutinos are transferred to the SD they are quickly locked into a 
mean motion resonance with Neptune (similar to the behavior of 
SDOs analyzed by Fernandez et al. d2"U04j l and Gallardo (2006). 
In the Centaur zone (this is the zone of q < 30 AU ) the distribu- 
tion of escaped plutinos is similar to that of SDOs in the Centaur 
zone obtained by Di Sisto & Brunini (2007). 

The orbital evolution of escaped plutinos in the Centaur zone 
can be grouped into the four dynamical classes proposed by Di 
Sisto & Brunini d20071 >. There are more particles that have the 
dynamical behavior of the second class and it is notable the great 
frequency of the presence of Kozai resonances in all the four 
classes. There are also several mean motion resonances densely 
populated in the ranges of 30 < a < 50 AU. 

The escaped plutinos have a mean lifetime in the Centaur 
zone of 108 Myr, greater than that of Centaurs from SD of 72 
Myr . Escaped-plutinos live more time than SDOs in the greater- 
perihelion Centaur zone, causing a slower diffusion to the inner 
Solar System of escaped-plutino orbits than of SDOs orbits. 

The present rate of injection of plutinos with radius greater 
than 1 km to the Centaur zone is between 1.6 to 16 plutinos ev- 
ery 100 years and the number of plutino-Centaurs with radius 
greater than 1 km would be between 1.8 X 10 6 - 1.8 x 10 7 . Both, 
the rate of injection and the number of Centaurs from plutinos 
are much less than the contribution from the SD obtained by 
Di Sisto & Brunini (12007ft . Then, plutinos would represent a 
secondary source of Centaurs and their contribution would be 
a fraction of less than 6% of the total Centaur population. 
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